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TECHNICAL MEMORANDUM ITO. 611 

MECHANICAL SIMILITUDE AND TURBULENCE* 
3y Th. v. Karrnan 

The development of hydrodynamics within the last decade 
has shown that skillful application of the equations from the 
dynamics of ideal fluids quite often brings clarity into such 
phenomena which in themselves are not independent of the vis- 
cosity. The vortex equations, in particular, proved themselves 
very useful. I may be allowed to mention the theory of the 
t vortex street by which we are able to reproduce the mechanism 

.* of. the form resistance with suitable approximation under stated 
, - H conditions, • although such a resistance is precluded in a fluid 

- • t . which is perfectly inviscid. Disregarding for the present the 
t 

origination of the vortex, the stream attitude in the wake of 
thes body may be described approximately correct by the repre- 
sentation of individual vortices, without transgressing the law 
governing the motion of such vortices in an ideal fluid. An- 
other striking, example is the theory of the induced drag of 
wings, which likewise shows the extent of applying the' vortex 
equations without overstepping the bounds of the dynamics of 
: ideal fluids. 

t 

* ^ But the prospects are ostensibly less promising for turbu- 

♦"ileclianische Ahnlichkeit und Turbulenz." Reprint fyorn Nach- 
richten von der Gesellschaft der Wissenschaf ten zu Gottingen, 
1930 - Fachgruppe 1 (Mathematik) No. 5, pp. 58-76. 
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lent fluid motion. It naturally is clear from the beginning 
that the balance between frictional and inertia forces prepon- 
derates in the problem of nascent turbulence-, so that the cited 
method does not enter into consideration at all in this case. 
But there remains the topic of the "fully developed turbulence," 
the problem according to the theory of "hydraulic stream atti- 
tude," which is perhaps of still greater importance for the 
practice. Still there are some indications which hold at least 
some promise of success in am attack on the problem. when disre- 
garding the frictional forces, or better expressed,, with the 
frictional forces confined to a definite zone, for instance, 
directly adjacent to the wall* Some time ago W. Fritsch* meas- 
ured in the Aachen Aerodynamic Institute the velocity distribu- 
tion in grooves between two parallel walls by constant wall dis- 
tance and for varying degrees of wall roughness. The results 
showed that the distribution curves - apart from the immediate 
neighborhood of the wall - are almost exactly superimposed as 
soon as the shearing stress at the wall assumes the same value, 
regardless of whether the fluid passes by a polished, smooth 
wall or, with correspondingly lower velocity, past a rough wall, 
or even a wall being saw-like in profile. The flow resistance 
in these grooves with very rough walls follows the so-called 
"square" law, i.e., it is proportional to the fluid density and 
to the square of the velocity, but unaffected' by the degree of 

fluid viscosity. Is it not feasible therefrom to surmise that 
*W. Fritsch, Zeitsch. f. angew. Math. u. Mech. , 1928, p. 215. 
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the viscosity is without predominating effect on the development 
of the velocity distribution? The present paper represents a 
method of attack in this sense, and endeavors an attempt to make 
the laws of turbulent flow in grooves amenable for calculation 
with a minimum of arbitrary assumptions. 

We are indebted to 0. Reynolds* for his explanation of the 
existence of turbulence as 'oscillatory motion, which - in con- 
trast to molecular unrest, responsible for the laminary fric- 
tion phenomena — he designated as molal fluctuating motion. 
Segregating the velocity components dependent on the time (We 
call them fluctuating, or oscillatory components) further the 
pressure variations from the mean values which correspond to 
the basic flow, the general hydrodynamic equations reveal the 
momentum components transmitted by the oscillatory motion as 
supplementary stresses from the standpoint of the basic motion, 
which are defined by mean values of products from the oscilla- 
tory velocities, 

I confine myself to the case of two-dimensional flow and 
express the mean of velocity x in the mean flow direction by 
U, that of velocity y by V, and the oscillatory components 
by u and v. Then the supplementary stress components read as 

° x = - P ? 

O y = ~ P - 

T = — P U~V 

*0. -Reynolds, Phil. Trans. A, Vol. 186, 1894, p. 123. ~ 
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These components are conformal to the so-called laminar 
stress components originating.. : in an unsteady, viscous fluid due 
to the molecular motion as given Toy 

"°y = 2 \i || (2) 

.. ' ' . t- =; \x m + m- 

The analogy becomes particularly . apparent when we include 
the derivation of the coefficient of viscosity. H, which is. r; 
defined up to a numerical factor, the product u = p X. c, 
where P = density, X = mean path length, and c = molecular 
velocity. . , 

Take, for instance, the case of. parallel, flow in direction 
x, that is, U = U(y) , V(y) = 0, and consider the derivation 
of the shear stress - first, conformal to 

t = - P u v 



and then, according to 



T =u *?s p x 0 dO , 

dy - . • dy 



As oscillatory components in the molecular motion, we have 

u = ± c ± I lH2. v = ± c\ and to the extent that 
dy 

u * ± o - I @, is coupled with,, v - c # .and u = ± c + I g 



c. Obviously,, uv = c 
analogy between the two expressions. 



with v = - Ct Obviously,, uv = c I which proves the 
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The next step to render the turbulent shear stress calcu- 
lable is to formally introduce a "turbulent friction coefficient, 
as proposed by Boussinesq* and notably, by Stanton.** The lat- 
ter, appreciating the velocity distribution developed from the 
center in a smooth pipe, disclosed the relation 




where t q = shearing stress transmitted to the wall, r = pipe 
diameter, and y = distance from middle of the pipe. This law 
of similitude is summarily adhered to by the introduction of an 
average turbulent friction coefficient M^, which is proportion- 
al to the quantity U max r.*** 

r = const. U max r |2 . 

H. F. Treer**** recently attempted to fit Stanton ! s theorem 
to the more up-to-date test material by having report to addi- 
tional empirical assumptions. It is, however, evident that the 
equation is carried too far in forming the mean value and is, 
on the face of it, inappropriate for interpreting actual veloc- 
ity distribution. Stanton's application may' be adjudged as a 
summary averaging examination of similitude; to conceive the 

mechanism of turbulent flo w atti tu de, the st udy of similitu de J 

*J. V. Boussinesq, Mem. pres. par div. .sav. , Paris 23, 1877; 

24, 1377. 

**T. E. Stanton, Proa. Roy. Soc, Vol. 85, 1911, p. 366. 
* ** /t 

To be more exact Mt = cons *- V tt r > since > however, - at 

least, for large index figures - -Jp is practically proportional 
*° U max ^, the above is likewise approximately correct. 
F. Treer, Phys. Zeitschr. , 1929, Vol. 30, p. 542. 
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must, to my way of thinking, be corded slightly different. 

Similitude of Oscillatory Motion ... 

Let us check the following assumption with: the mean flow 
to be a parallel flow in direction x; 'the oscillatory flow to 
consist in disturbances, of relatively, limited extent in direc- 
tion y and - at least, within a certain time interval - to be 
carried along -by the main flow approximately as steady flow con- 
figurations.- Our problem' shall be to find under what condi- 
tions these flows can become similar among each other, so that 
the flow attitudes neighboring two point's,' which conform to 
different y values, vary only by a multiplicative factor of 
the oscillatory velocity and, in the length* measure of the field 
of flow. In other words, we assume similitude , pf the oscilla - 
tory attitude, irrespective of the location of the point in 
wh ose vicinity the oscillation is examined . 

We so place the coordinates that this pertinent point .falls 
on axis y = 0, and develops the mean velocity according to y: 

U - U«p7 ; + U» 0 . .... . " 

Then we write the stream function adjacent to' the point 

* (x,y) = U< 0 % + U» 0 £ ■+ ... + ♦ (^y), (5) 

where * (x,y<) becomes tire stream function of the oscillatory 
motion. 

Now it is desired that only I and A be affected by the 
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point, that is, by U« 0 , U" 0 , etc, but that function f(t,Ti) 
he unaffected when 

x = i e 

y. • = l r\ (6) 

f ' = a f (e,ri) 

The hydrodynamic equations for plane steady flow may be 
combined in the so-called vortex transfer equation 

ai _ at aM = v> a -a $ (7) 
ay sx ax ay 

where 

3.x 2 • 9y 3 

is the vortex intensity, and v the kinematic viscosity. Dis- 
regarding the friction terms and limiting ourselves, in accord 
with the assumed definition of the field of oscillation in the 
vicinity of axis y = 0. : to the first digits .-in ;the terms in- 
duced by the principal, motion, we obtain ■ 

Cut. y + 2T) m - |l u«o - 2* ^ = o (8) 
Introducing (6), we have 

TT , , a 3 Af a af TT „ . a s /af aAf ?tf aAf> _ 0 fg x 

where symbol A now pertains to the variables i and r\ . This 
equation is resolvable independently of A, I and U' 0 , u""o, 
providing 
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UV I ~| 



U" 0 = — 



(10) 



or 



and 



u» 



O. 



* s ^4 ~ i * «'p . (id 

u O 

To denote the significance of these relations, we now com- 
pute the shearing stress t: 



or 

t = P I 2 U'o 2 (13) 

The conditions of similitude may then be combined as follows: 
a) The field of oscillation retains a length indicative 
of the length measure of the disturbances, defined by 

• • dU, ■ 

U'o-..: <iy : .. . ■ : 



U»o d!_u 

dy 2 ' • 

b) The shearing force is proportional to the density, to 
the square of the characteristic length i and to the square of 
the velocity gradient 

'. dy v.r.: 

Having defined length I to only one factor, we write 

r = p i* (14) 
Vdy/ KJ -*' 
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and I - k 51 (15) 

where k is a nondimensional constant, dependent only upon the 
nature of the fluctuating mechanism, and upon the solution of 
the equation for f(£,T)). It is the sole constant entering the 
here developed theory of turbulence, 

Prandtl's Theory of "Mixing Length" 

Our data may for the present serve as corroboration of the 
suggestive equation put forth by Professor L/'Prandtl* to demon- 
strate the laws of turbulent flow. Ke proceeded with the for- 
mula for the shearing stress 



T = - P U V 



and evaluated u = I analogous of the equations of the molec- 

dy 

ular theory, by introducing the length I as the nixing length 
or distance (analogous to the idea of the mean free path in the 
kinetic theory) . On the premises that a particle in the oscil- 
latory motion travels the distance I perpendicular to the basic 

flow without momentum interchange p I § becomes, in fact, the 

dy 

momentum which the particle transfers to; the layer into which it 
Y^as displaced by the mixing motion, Prandtl then assumed the 
fluctuating component v proportional to the mixing length i 
and to the absolute value ~\ of the velocity gradient, so 
that in agreement with (13), t = p l 2 

Viewing the matter from the viewpoint of the oscillatio n 

♦Compare, for example, page G2 of "Proceedings of the Second In- 
ternational Congress for Teolinisohe llechanilc, 1937." 
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theory, this assumption signifies that the correlation between 
u and v is supposedly unaffected "by the location. This hypoth- 
esis is so obviously like our assumption of similarity, that -it 
is not surprising to find the equation confirmed. 

But now our study of similitude reveals decidedly more, for 
it yields an equation defining the mixing distance which in 
first approximation appears as ratio of the first two differen- 
tial quotients of the basic velocity. 



Consider the basic flow in direction x between two paral- 
lel walls y = ± h (Fig. l). With t 0 as shearing stress at 
the wall, it equals 

at y distance away from the channel center, so that the veloc- 
ity distribution becomes 



Flow between Parallel Walls 




(17) 



or 





= k 



One integration yields 




or (a = a constant) 
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u * = h-r^ — ^7=- (18) 

2k Jh a - JT 

The constant a, is to be defined from the relation with 
the conditions at the wall. For great Reynolds Numbers ^ as- 
sumes a high value near the wall and finally approaches the lim- 
iting value = Ifi which, due to the smallness of \x re- 
dy ^ 

ferred to ~ at a distance away from the wall, is very high. 

rl T 7 

As a result the value y, at which ^ becomes infinite, may 
be allowed to coincide with y = h, so that 

1 V 0 



dU 



wherefrom integration yields the velocity distribution. The 
maximum U max is reached at y = 0. Then 

max o V P 2kh fy 



(19) 



or 



u = u max + [log (i - /T) + /T 

This formula was compared with two measurements of Donch and 



(20) 
(20a) 



U 



Nikuradse.* Both defined the values -p~, leaving k as sole 
constant to be determined, which appears to have a value of 

about k = 0.36. The resistance law yields, as shown l ater, 

*F. Donch, Forschungsarbeiten, No. 282, 1926. 
J. Nikuradse, Forschungsarbeiten, No. 289, 1929. 
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k = 0.33. j 

Figure 2 illustrates the comparison, which is found to be 
good. One more word about the* behavior of mixing length I . 
The calculation yields • '* 

i = k g=3 k h/F(!-/0- ..; (21) 

Near to the wall we may write , ...... 

SO that , v 

i - 2k h {/i44r^-w,(i I k^- •••)• < 22 > 

The course of ^ = f , is noted in Figure 3. It does not 
correspond to Donch and -Kikuradse 1 s data -(dotted" line) , which 
they arrived at with formula • • -.-w;* • r ■■ 



dU 



Their chief contention is 5 , that ; quantity I assumes a constant 

value in the channel center, whereas it attains a maximum at 

■ ■ ■ 

y i " f 11 and - dro P s *° zer 9 ^or ; y x =.h; according to our calcu- 
lation. Cur comments are: 

In the first place, '1 is exceedingly difficult to deter- 
mine on account of the uncertainty m of the differentiation of 
the point-by-point recorded velocity curve; then they compute 
the mixing length i n . the.: channe 1 c enter - by - a mo dlf i ed formula, 
the authenticity of which is not" quite "beyond question. It is 
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assumed forthwith that $1=0 for y = 0, that is, the veloc- 

dy 

ity is smoothed out at the apex. Oixr formula shows a hump for 
y = 0, a sure sign that the degree of approximation in the cal- 
culation is insufficient at this point, but it is strange that 
the unsmoothened velocity curve also shows a definite bend* 
I believe the observation material was not extensive enough to 
define the behavior of I in the channel center conclusively. 

The Resistance Law 

By a high velocity gradient f|2, as near to the wall, our 
similitude consideration no longer holds true, because the omis- 
sion of the viscosity appears to be no longer justified. It is 

a question of whether \x ^ can be disregarded or not along 

dy 

side of the kinematic shearing stress - P u v. It has been 

proved experimentally that near to the wall something like a 

"laminar layer" exists and, in order to define the resistance 

law, i.e., the U^g^ value corresponding to a certain value of 

-fi, the connection from the velocity curve to the laminar layer 
P 

must be established. 

Before discussing a more exact theory of the oscillatory 
field, we attempt two different applications which, however, 
yield the same plausible results; 

a) We assume the mixing length I to diminish below zero 
to a value proportional to the thickness of. the laminar layer. 
This minimum is supposedly reached at the boundary of the lami- 
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nar layer where the lardnar layer is then attached. 

To carry this idea, through we resort to G. J. Taylor^ con- 
ception, who arrived. at the. conclusion thq,t the vortex division 
(i.e., the. mixing length I.) of the tiirbulenoe. at the wall can 
only be affected "by the ; shearing stress t q? and that it must 

r -- j ------ - lor 



be proportional., to the value 



Writing I = a 



"ITi. 



this 'limiting value, 'we "obtain for the relevant " quantity of- 
1 ~ P& r-e 



and U = - U^ ax U(y.) , that is for the velocity 



y = i ay 
' k > 0 

difference "between channel center and laminar layer "boundary 



tje _ 1 / V 



log 



•+ const . 



(23) 



T ee 



Thp.n we estimate the difference -in velocity U between 
both boundaries of the laminar layer , • i. e. , between the wall 

and the free boundary. Since the thickness on one side is 

T ee - • ' " '-' •'• v 



6 2 



v 



<. ...and we can put t c = \x - 



IT 



on the other, we obtain 



~ ' and the difference between the wall and the .channel 



center beeves 



Umax = U e + U ee = £ 



T 



log 



v 



(24) 



The . oonstent^ ; ; -^nd A ■; are ; independent .of dimensions and Roy- 
nolds ITumber . » ... 

\/tr) .'The^result is identical when, we assume the "turbulent 



coiitpQn.Snt".;.'T - -jjt., 



au 



to. be • proportional - to P. T 



in- 



ay - : \ dyy 

«tead of i the whole;, shearing.- stress.. Then I r.aay„ quietly dimin- 
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ish to zero and the velocity distribution near to the wall be- 
comes 



This is the same equation used by Ur. Wada* in formulating 
the resistance law, but which, owing to having been published 
in an inconspicuous place, has not received sufficient attention. 
Unfortunately, Mr. Wada held the equation valid for the whole 
channel which made his formulas a little too complicated, al- 
though the results in this paragraph are really contained in 
his report in an implied form. 

To der.ive the resistance law, I slightly deviate from the 
conventional parameters and introduce the factors 

■p u max k 
K m 77— 



* = / (26) 
p *nax 



In place of the mean velocity the maximum velocity appears 
as reference quantity, so that (34) yields 



log (R m JH) + A - i log 2 

(27) 

or 

= log (R m /T) + c" (28) 

We -derived the relation "between U max and ^ for plane flow. 
♦Journal o^ _ Japane"sV^oc". Naval~Arch. 41, 1927, p. 103. ~ 
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But it can also be proved applicable to circular pipes, where 
even constant k retains the same value and where constant C 
alone is different. In our derivation we used the mixing dis- 
tance k x wall distance and estimated the thickness of the 
laminar layer. But it may be assumed that this is valid for a 
circular pipe also, for the further course of the mixing dis- 
tance conies in evidence only in the 0 constant* 

Figure 4 shows some observations of ^=^= and log 10 R m 
U r 

(R m = , gjp - ? r = pipe 'diameter).* It will be noted that there 
is a linear relation between both quantities, extending from 
R m = 3000 to R m - 1,600,000. Constant 

k.= 0*38 

(29) 

0 - 1.83 

the first being perhaps universal, the second applying particu- 
larly to circular pipes. 

The So-called Power Laws 

It is known that the resistance law within large ranges 
of Reynolds Numbers can be adequately expressed by the interpo- 
lation formula 

, a const. (30) 
R m 

Upon this premise it can be proved by a line of reasoning 
advanced by Dr. Prandtl that the velocity distribution (figured 
as beginning at the wall) is given by formula 



♦The data on very large Reynolds Numbers were supplied by Mr. 
Hikuradse, who placed his, as yet unpublished material, at my 
disposal.- 
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Ho * 

u = const. J 2$ (ijLl) 



(31) 



The validity of this power formula ceases in the immediate neigh- 
borhood of the wall (effect of laminar layer), "but extends in 
surprising manner to the other side nearly to the channel cen- 
ter. The exponent h drops as the Reynolds Number increases 
(within about 1 to in the range examined thus far) . 

This puzzle is simple to explain. Slightly transformed, our 
derived formulas express the 



Resistance law: 

Velocity distribution: 

So when we make n = 
written as 

Resistance law: 
Velocity distribution: 



-L= = a + b log (RjT) 

a/tjTT ... 



J<2 



(32) 



h= - a ' + b ' log 



m 



3 - m 



, equations (30) and (31) can be 



= const. (Rj~V) n 



U = const. (^_L) n 



(33) 



p 



The decrease in exponent n with the Reynolds Number becomes 
readily apparent from the following: 

Comparing flow attitudes in the same channel which corre- 
spond to an identical value of t Q9 and permitting the viscos- 
ity to vary, for instance, v> to decrease gradually, the veloe- 
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ity distribution is obviously the same up to the point where 
the influence; of^ta There a kind 

of laminar -layer is -set fug*.- although. ^ .. . 

shift along ; the v e lo c i t y ; di s t r ibu t i o n, curve as.. the viscosity, de- 
creases, i.i e-. y as ;thejReynoldsf Ku^iber .rises. ^The . conjoint ed. 
distribution.: curve fr-, .-. ; . ; . v n > I ... ... 

which is, aside from the channel center.^ -exponential . (Fig. 5) 5 
is to be approximated , to ..y U 1/n ;■ by a series of power curves 
which touch the axis y = 0, inducing the contact point to grad- 
ually shift -toward the right; accordingly l/n increases and 
n decreases. \ 

< f .,<\ ■ \ ..... 0 u g h n e s s 

It has been repeatedly emphasized that the ratio ; between 
the thickness of 'the laminar^ layer . and the mean protuberance 
of the roughness predominates the phenomena on rough walls. I 
call a wall Tough when . the ;pro j ections are large referred to. the 
laminar layer thickness. In this case it nay be assumed that 
the • mininuni value of the.; mixing, -length I is no t. .conditioned by 
the thickness, but by the size^.of :the roughness. ..elements.. Thus 
if - •€ -is- -the- mean roughness projection (the characteristic 
length measure of the rQughne&s) y .the. .minimum yalue, of I may be 
made proportional to c r which discloses as relation between 
velocity in channel center U max and the shearing force at the 
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wall 



u max 



i /^"(log | + const.) (34) 



or for the resistance coefficient 

k a/~2~ h 

/rn— 1 = lOg - + Const. 

In other words, the flow resistance follows the "square" 
law, and the resistance coefficient is dependent on the rela- 
tive roughness according to (35). This equation admits of a 
check when comparing experiments in grooves, where the distance 
of the walls (2h) varies, but where the nature of the walls 
(quantity c) is to remain constant. This, however, is to be 
treated in a future report. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 3,4,5 
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